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Abstract Magnetars are proposed to be peculiar neutron stars which could power their X-ray 
radiation by super-strong magnetic fields as high as > 10^^ G. However, no direct evidence 
for such strong fields is obtained till now, and the recent discovery of low magnetic field 
magnetars even indicates that some more efficient radiation mechanism than magnetic dipole 
radiation should be included. In this paper, quantum vacuum friction (QVF) is suggested to 
be a direct consequence of super-strong surface fields, therefore the magnetar model could 
then be tested further through the QVF braking. Pulsars’ high surface magnetic field inter¬ 
acting with the quantum vacuum result in a significantly high spindown rate ( P ). It is 
found that QVF dominates the energy loss of pulsars when pulsar’s rotation period and its 
first derivative satisfy the relationship ■ P > 0.63 x where ^ is the ratio of 

the surface magnetic field over diploe magnetic field. In the “QVF -f magnetodipole” joint 
braking scenario, the spindown behavior of magnetars should be quite different from that in 
the pure magnetodipole model. We are expecting these results could be tested by magnetar 
candidates, especially the low magnetic field ones, in the future. 
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1 INTRODUCTION 


Kinematic rotation was generally thought to be the only energy source for pulsar emission soon after the 
discovery of radio pulsars until the discovery of accretion-powered pulsars in X-ray binaries. However, 
anomalous X-ray pulsars/soft gamma-ray repeaters (AXPs/SGRs, magnetar candidates) have long spin pe¬ 
riods (thus low spindown power) and no binary companions, which rules out spin and accretion in binary 
system as the power sources. The first SGR-giant flare was even observed in 1979 jMazets et al. 1979b . and 
Paczynski (1992) then pointed out that the super-strongmagnetic held may explain the super-Eddington 
luminosity. AXPs and SGRs are thereafter supposed to be magnetars, peculiar neutron stars with sur¬ 
face/multipole magnetic fields (10^''^ G ^ 10^® G) as the energy source, while the initially proposed strong 


dipole fields could not be necessary (e.g., Tong et al. 20131. Moreover, the discovery of low magnetic held 
magnetars (IZhou & Chen 201^ IRea et al. 20101 IRea et al. 20121 IScholz et al. 2012t in recent years indi- 
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Besides failed predictions and challenges in the magnetar model (IXu 20071 Tong & Xu 201 1| |, one of the 
key points is: can one obtain direct evidence of the surface strong fields? Here we are suggesting quan¬ 
tum vacuum friction (QVF) as a direct consequence of the surface fields, and calculating the spindown of 
magnetar candidates with the inclustion of the QVF effect. 

Magnetodipole radiation could dominate the kinematic energy loss of isolate pulsars (e.g., 


Manchester & Taylor 1977|IDai & Lu 19981 Lyubarsky et al. 2001 IMorozova & Ahmedov 20()8] |. The de¬ 

rived braking index n = (H is the angular velocity of rotation) of a pulsar is expected to be 

3 for pure magnetodipole radiation. As a result of observational difficulties, only braking indices n of a 
few rotation-powered pulsars are obtained with some certainty (lYue et al. 20071 [Livingstone et al. 2007| 
http://www.atnf.csiro.au/research/pulsar/psrcat/). They are PSR J1846-0258 (n = 2.65 ± 0.01), PSR 
B1509-58 (n = 2.839 ± 0.001), PSR Ill 19-6127 (n = 2.91 ± 0.05), PSR B0531h- 21 (the Crab pulsar, 
n = 2.51 ±0.01)), PSR B0540-69 (n = 2.140±0.009) and PSR B0833-45 (the Vela pulsar, n = 1.4±0.2). 
These observed breaking indices are all remarkably smaller than the value of n = 3, which may sug¬ 
gest that other spin-down torques do work besides the energy loss via dipole radiation (jXu & Qiao 200 Ij 


IBeskin et al. 1984llAhmedov et al. 2012IIMenou & Perna 20011 Contopoulos et al. 2006 Alpar et al. 2001 


IChen et al. 20061 IRuderman 20051lAlIen et al. 1997llLin et al. 20041 [Tong & Xu~20T4l|Tong 2015^ 

Recently, the research of Davies et al. shows that the QVF effect could be a basic electromag¬ 
netic phenomenon (IDavies 20051 ILambrecht et al. 19961 jPendry 1997| [Feigel 2004| [Tiggelen et al. 2006| 


Manjavacas et al. 20101. If the quantum vacuum friction exists, the dissipative energy by QVF would cer¬ 


tainly be from rotational kinetic energy of pulsar. The loss of rotational kinetic energy of pulsar by QVF 
may also transform into pulsar’s thermal energy or the energy of pulsar’s radiating photons which might 
not be isotropic. This is the same argument as in the work of Manjavacas et al. (2010), in which the au¬ 
thors argue that at zero temperature, the friction produced on rotating neutral particles by interaction with 
the vacuum electromagnetic fields transforms mechanical energy into light emission and produces parti¬ 
cle heating. Pulsar may transfer its angular momentum to the vacuum when pulsars rub against quantum 
vacuum since the angular momentum is conserved. In this case, vacuum may work as an standard medium 
( Dupays et al. 2008| . Dupays et al. (2008, 2012) even calculated the energy loss due to pulsars’ interaction 
with the quantum vacuum by taking account of quantum electrodynamics (QED) effect in high magnetic 
field. The calculations indicate that when the pulsars’ magnetic field is high, QVF would also play an im¬ 
portant role to cause the rotation energy loss of pulsars. Thus, it is necessary to take QVF into the rotation 
energy loss of pulsars, especially for highly magnetized pulsars on surface, like magnetars. 


In this paper we assume that pulsar interacts with quantum vacuum as in the work of Dupays et al. 
(2008) and consider the difference between the surface/toroidal magnetic field and dipole/poloidal magnetic 
field. The braking indices for pure QVF radiation and surface magnetic field of magnetars for the “QVF ± 
magnetodipole” joint braking model are calculated. 

The paper is organized as following. After an introduction, we deduce the relation between the dipole 
magnetic field and the braking index of magnetars in the second section. The calculated results and analysis 
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2 SPINDOWN AND BRAKING INDEX OF MAGNETARS 


A pulsar has the power of magnetodipole radiation of 

2 


where 


-Edip = -3C 


1 ^ 

n = -Bdip-R Sin 0 


( 1 ) 


( 2 ) 


is magnetic dipolar moment and c is the speed of light in vacuum, Bdtp is the dipole magnetic field, R is the 
pulsars’ radius, 9 is the inclination angle. For general pulsars, surface magnetic field approximately equal 
to dipole magnetic field because multipole magnetic field attenuate to little. However, for magnetars there 
is a surplus of attenuate multipole magnetic field as its extraordinarily strong surface magnetic field. So 
magnetars’ surface magnetic field Bgurf include dipole magnetic field B^ip and multipole magnetic field. 
We suppose that the ratio of surface magnetic field and dipole magnetic field 

.^surf 


e = 


^dip 


(3) 


is a constant. The pulsar rubs against the quantum vacuum and then loses its rotation kinetic energy 
( |Dupays et al. 2008] ! of 


.^qvf — tr - 


Stt sin^ 9 Bi 


surf 


(4) 


16 cBl P2 ’ 

where a = jhe ~ 1/137 is the coupling constant of electromagnetic interaction. Be = 4.4 x lO^^G is 

the QED critical field and P = 2Tr/His the spin period. 

The pulsars’ typical radius R = 10®cm is adopted. Set inclination angle 9 = 90° for the sake of 
simplicity. Considering the relation (|2l) between the magnetic moment of pulsars and magnetic field in polar 
region of pulsars, we can obtain the ratio of the energy loss due to QVF over that due to magnetodipole 
radiation 

Eqv! 


E, 


dip 


= 7.69 X 


(5) 


Assuming the pulsars’ rotation energy loss coming from both magnetodipole radiation and QVF, i.e. 
E = Edip + Eqv{, the total energy loss of pulsars are given by 

Stt sin^ 9 


Ec^-- 


-e- 


( 6 ) 


3c3 16 cB2 p2 

From the pulsars’ rotation energy loss E = IHH, where I is the inertia of momentum with typical value 
/ = lO'^^g cm^, we can obtain a relationship between pulsar’s period and the period derivative with respect 
to time 


P = 


27r^i?® sin^ 03Qfi?^ sin^ 0 


3c3/ P GAttBc^Ic 
Using the relation of H and P, the braking index can be obtained 

1 / 27r^i?® sin^ 0 T^dip 3a i?^sin^0 


BipP^- 


Ic^ 


6471 B/Ic 


Bt,pPe 


(7) 


( 8 ) 


Numerically, the braking index can be written as 
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where 


(10) 


with i?dip,i 2 = 10 ^^ildip- We can also express the ratio of the energy loss due to QVF over that due to 
magnetodipole radiation by pulsar’s period(P) and period derivative(P) from equation (|5]l and (|7]i 


= 7.69 X 10-24 




dip 


8768 


8768 




8768 


3P4 


( 11 ) 


Numerically, the above equation can be written as 


:^ =-i + wi+ 3.16 X 10i6^4p3p_ (12) 

F^dip 2 V 4 

3 THE NUMERICAL RESULTS 


The periods of observed pulsars are distributed mainly in the range from 0.1s to 5s (The ATNF Pulsar 
Catalogue: http://www.atnf.csiro.au/research/pulsar/psrcat/). Using Eq. 0 we plot the ratio of Pqvf/i^dip, 
as a function of the period P in Fig. 1 for ^ = 10 and in Fig. 2 for ^ = 100. From Fig.l we can see that 
QVF may play an important role when the dipole magnetic field is higher than ^ lO^^G for pulsars whose 
period are between 0.1s and Is. Most of observed pulsar’s magnetic field derived from pure magnetodipole 
radiation are in the region lO^^ — lO^^G, however, if QVF is included in pulsars’ energy loss, the derived 
magnetic field could be lower. Thus it is necessary to independently measure the magnetic field of pulsars 
so that we can judge whether QVF has important contribution to pulsars’ rotation energy loss. 

From Fig.2 we can see that QVF may play an important role when pulsars dipole magnetic field i^dip > 
104° for most pulsars’ braking. For millisecond pulsars the derived magnetic field from magnetodipole 
radiation is already so low {B^ip < lO^^G) that we can neglect the QVF’s contribution to its rotation 
energy loss, but for magnetars the derived magnetic field from QVF is already so high (Bdip > IQi^G) 
that we have to consider the QVF’s contribution. 'We can also express the ratio of the energy loss due to 
QVF over that due to magnetodipole radiation by pulsar’s period(P) and period derivative(P) as shown in 
Eq. m- Erom this equation we can obtain that QVF dominates the energy loss of pulsars when pulsar’s 
rotation period and its first derivative satisfy the relationship P^ • P > 0.63 x where ^ is the 

ratio of the surface magnetic field over diploe magnetic field. According to above relationship and current 
observed data for confirmed magnetars (see Table 1) QVF will dominate the rotation energy loss in all of 
the magnetars’ spindown. 

Substituting the observed value of P and P into Eq. 0, the magnetic field of pulsars can be calculated. 
We compute the currently confirmed magnetars’ magnetic field and list the results in the last column PJj"p 
of Table 1. The fourth column P^ip is derived from pure magnetodipole radiation. The calculated results 
manifest that the derived dipole magnetic field Pdip from pure magnetodipole radiation is about 10^(^ = 
10) and 104(^ = 100) times larger than PJj"p obtained by combining QVE and magnetodipole radiation. 
And the derived surface magnetic field Psurf from pure magnetodipole radiation is about 100 times larger 
than P)j"p inferred by combining QVE and magnetodipole radiation for both ^ = 10 and ^ = 100. 

If P = pQVF, from Eq. @ we can obtain 17 cx therefore braking index n = 1 for pulsar’s spindown 
by pure QVF. Eq. 0 show that pulsar’s braking index is between 1 ^ 3 in the ‘QVE + magnetodipole” 
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P(S) 


Fig. 1 The ratio of a pulsar’s energy loss rate from QVF over that from magnetodipole radiation, 
as a function of period, where ^ = 10. 



P(s) 


Fig. 2 The ratio of a pulsar’s energy loss rate from QVF over that from magnetodipole radiation, 
as a function of period, where ^ = 100. 

so magnetars have bigger /(Sdip, PA) function value (see Eq. ( fTOb ) which result in QVF dominating 
magnetars’ braking and its braking indices being about 3. However, for some low magnetic field millisecod 
pulsar, minor f{Bdip, PA) function value lead to magnetodipole radiation becoming main energy loss 
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Table 1 The parameters and the inferred magnetic field of magnetars. The magnetars’ data 
of period (P), the period derivative (P) and dipole magnetic field (Sdip) are from the McGill 
SGR/AXP Online Catalog (http://www.physics.mcgill.ca/'^pulsar/magnetar/main.html). The last 
column of table, is inferred magnetic field from our model based on both magnetodipole 
radiation and QVF. 


Name 

P(s) 

P(10"“s/s) 

Pdip(lOi^G) 

Pi-?f(10“G, C = 10) 

e = 100) 

CXOU JO 10043.1-721134 

8.020392(9) 

1.88(8) 

3.9 

5.946 

5.946 

4U 0142-t61 

8.68832877(2) 

0.20332(7) 

1.3 

3.342 

3.342 

SGR 0418-1-5729 

9.07838827 

< 0.0006 

< 0.075 

< 0.770 

< 0.770 

SGR 0501-1-4516 

5.76209653 

0.582(3) 

1.9 

4.817 

4.817 

SGR 0526-66 

8.0544(2) 

3.8(1) 

5.6 

7.082 

7.082 

IE 1048.1-5937 

6.4578754(25) 

2.25 

3.9 

6.565 

6.565 

IE 1547.0-5408 

2.06983302(4) 

2.318(5) 

2.2 

8.791 

8.791 

PSR11622-4950 

4.3261(1) 

1.7(1) 

2.7 

6.766 

6.766 

SGR 1627-41 

2.594578(6) 

1.9(4) 

2.2 

7.905 

7.905 

CXO J164710.2-455216 

10.6106563(1) 

0.083(2) 

0.95 

2.541 

2.541 

IRXS J170849.0-400910 

11.003027(1) 

1.91(4) 

4.6 

5.516 

5.516 

CXOU 1171405.7-381031 

3.825352(4) 

6.40(5) 

5.0 

9.718 

9.718 

SGR 11745-2900 

3.76363824(13) 

1.385(15) 

2.3 

6.655 

6.655 

SGR 1806-20 

7.6022(7) 

75(4) 

24 

15.145 

15.145 

XTE 11810-197 

5.5403537(2) 

0.777(3) 

2.1 

5.229 

5.229 

Swift 11822.3-1606 

8.43772106(6) 

0.00214(21) 

0.14 

1.078 

1.078 

SGR 1833-0832 

7.5654091(8) 

0.439(43) 

1.8 

4.194 

4.194 

Swift 11834.9-0846 

2.4823018(1) 

0.796(12) 

1.4 

6.430 

6.430 

IE 1841-045 

11.7828977(10) 

3.93(1) 

6.9 

6.494 

6.494 

3XMM J185246.6-I-003317 

11.55871346(6) 

< 0.014 

< 0.41 

< 1.594 

< 1.594 

SGR 1900-1-14 

5.19987(7) 

9.2(4) 

7.0 

9.856 

9.856 

IE 2259-1-586 

6.9789484460(39) 

0.048430(8) 

0.59 

2.466 

2.466 

PSR 11846-0258 

0.32657128834(4) 

0.7107450(2) 

0.49 

10.379 

10.379 


magnetodipole radiation, we use Eq. (|9]) to calculate the braking indices of magnetars. The results show 
that all the magnetars’ braking indices are around 1 for both ^ = 10 and ^ = 100. In the future, the model 
could be tested by comparing the calculated results to observed braking indices. This comparison can also 
provide further information to understand QVF. 

4 CONCLUSIONS AND DISCUSSIONS 

We investigate pulsar’s rotation energy loss from QVF and compare it with that from magnetodipole ra¬ 
diation in the different magnetic field range and different period range. We find that if the ratio ^ of the 
surface magnetic field over dipole magnetic field is fixed to 10(100), QVF could play a critical role for 
pulsars’ braking when Bsurt ■ P > 10^^(10^°)G-s, while it can be ignored when Ssmf ■ P < 10^°(10®)G-s. 
Magnetars may have high surface magnetic field and long period (Bsurf ■ P ^ lO^^G s) if the value of 
magnetic field is inferred by pure classical magnetodipole radiation. Therefore it is necessary to consider 
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We consider the difference between the surface magnetic field and dipole magnetic field of pulsars and 
compare the energy loss rate of pulsars due to magnetodipole radiation to that due to QVF. The results show 
that when a pulsar has a strong magnetic field or a long period (Bsurf ■ P > lO^^G-s for ^ = 10, i?surt ■ P > 
10^°G s for ^ = 100), comparing to QVF, the energy loss by magnetodipole radiation can be ignored, while 
when pulsars have weak magnetic field or short period (ilsurf • P < 10^°G-s for ^ = 10, i?surf • P < 10®G-s 
for ^ = 10) the QVF can be negligible. We consider that rotation energy loss of magnetars is the sum of 
the energy loss due to QVF and that due to magnetodipole radiation. Based on this joint mechanism of 
energy loss, the surface magnetic field of magnetars and braking indices are calculated. Our work indicates 
that when QVF is included in the process of rotation energy loss, the surface magnetic field of magnetars 
is 10 — 100 times lower than that in pure magnetodipole radiation model. In this joint braking model 
QVF dominates the energy loss of pulsars when pulsar’s rotation period and its first derivative satisfy the 
relationship ■ P > 0.63 x 10“^®^“^s^, where ^ is the ratio of the surface magnetic field over diploe 
magnetic field. Also, we obtain the braking index of magenetars is around 1 in the joint braking model. 
The efficiency of rotation energy losses generated by QVF in magnetars is very high compared to magnetic 
dipole radiation. Smaller magnetic field can generate a greater rotation energy loss by QVF comparing to 
magnetic dipole radiation. This may explain why magnetars which have great X-ray luminosity and low 
magnetic field (IZhou & Chen 2014llRea et al. 20101 IRea et al. 2012IIScholz et al. 20121 . 

We are expecting the results presented could be tested by X-ray observations of magnetar candidates, 
especially for the low magnetic field ones. X-ray data accumulated in space advanced facilities could show 
both timing and luminosity features for magnetars, and a data-based research would be necessary and in¬ 
teresting. Summarily, further observations for magnetars in the future would test our joint braking model as 
well as help us understand QVF in reality. 
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